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Abstract. Fortran subroutines to calculate helicity amplitudes with massive spin-3/2 particles, such as 
massive gravitinos, which couple to the standard model and supersymmetric particles via the supercurrent, 
are added to the HELAS (HELicity Amplitude Subroutines) library. They are coded in such a way that 
arbitrary amplitudes with external gravitinos can be generated automatically by MadGraph, after slight 
modifications. All the codes have been tested carefully by making use of the gauge invariance of the 
helicity amplitudes. 
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1 Introduction 



Gravitinos are spin-3/2 superpartners of gravitons in local 
supersymmetric extensions to the Standard Model (SM). 
If supersymmetry (SUSY) breaks spontaneously, graviti- 
nos absorb massless spin- 1/2 goldstinos and become mas- 
sive by the super-Higgs mechanism. Therefore, the grav- 
itino mass is related to the scale of SUSY breaking as well 
as the Planck scale like 



m 3/2 



(M SU sy)7^pi- 



(1) 



This implies that the gravitino can take a wide range of 
mass, depending on the SUSY breaking scale, from eV up 
to scales beyond TeV, and provide rich phenomenology in 
particle physics as well as in cosmology pQ. 

Although the gravitino can play an important role even 
in collider signatures when it is the lightest supersymmet- 
ric particle (LSP), there is few Monte Carlo event gen- 
erators which can treat themQ In this paper, we present 
new HELAS subroutines [3 for the massive gravitinos and 
their interactions based on the effective Lagrangian below, 
and implement them into MadGraph/MadEvent (MG/ME) 
v4 [H(5J|B] so that arbitrary amplitudes with external grav- 
itinos can be generated automatically^ 
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1 The SM with 
WHIZARD [2] . 

2 The Fortran code for simulations of the massive gravitinos 
is available at the KEK HELAS/MadGraph/MadEvent Home 
Page, [http : /7madgraph . kek . jp/KEK/] 



The effective interaction Lagrangian relevant to the 
gravitino phenomenology is [HEH] 

i 



\/2M P | 

x [trf-fPLf (dm* - pp R rY^ (dm 

- ^-f-fPiif (dm* + pPL-f-f^ (dm] 
-'^,Y]^X^ a F^\ (2) 

BMpi 

where ip^ is the spin-3/2 gravitino field, f % and (f> 1 are 
spinor and scalar fields in the same chiral supermultiplet, 
Pr/l = 3(1 i 75) is the chiral-projection operator, and 
Mpi ee Mpi/V&x ~ 2.4 x 10 18 GeV is the reduced Planck 
mass. The covariant derivative is 

D,=d^ + ig.TfA^ + igT^W* + ig'YB^ , (3) 

where g s , g and g' are the SU(3)c, SU(2)l and U(l)y 
gauge couplings, respectively, and Tg, TJf and Y are the 
generators of the SU(3) C (a = l,--- ,8), SU(2) L (a = 
1,2,3) and U(l)y groups. The field-strength tensors for 
each gauge group are 



p{3)a 



9 s j 3 
W - gff 
F^ a = d„B„-d„B 



abc j^b j^c 



i^( 2 ) a = d w a 



(4) 
(5) 
(6) 



and the corresponding gauginos A( a=3 ' 2,1 ) a are gluinos 
(g a ), winos (W a ) and bino (B), respectively. 

The paper is organized as follows: In Sect. [2] we give 
sample numerical results. Sect. [3] presents our brief sum- 
mary. In App. [21 we give the new HELAS subroutines for 
spin-3/2 particles, and in App. [B]wc describe how to im- 
plement the amplitudes into MG. 
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2 Sample results 

In this section, we present some sample numerical results, 
using the new HELAS subroutines, which are presented in 
Appendix [5] and the modified MG, which is described in 
Appendix [Bj 

In the gauge mediated SUSY breaking scenarios, the 
gravitino is often the LSP, and its phenomenology depends 
on what is the next-to-lightest supersymmetric particle 
(NLSP). Here we consider the stau NLSP scenario as well 
as the ncutralino NLSP one. 



2.1 Stau NLSP 

As a sample result for the stau NLSP scenario, we consider 
radiative f decays, 

f-^r-G 7 . (7a) 

Here we regard the stau as a purely right-handed stau 
for simplicity. Feynman diagrams shown in Fig. [T] and the 
corresponding helicity amplitudes are generated automat- 
ically by the modified MG. To study the spin-3/2 nature of 
the gravitino, we compare the G LSP case (|7a|) with the 
X? LSP case, 

Tr^t-x i1, (7b) 

where only two decay diagrams contribute; see Fig. [5J 

We evaluate the amplitudes for the both cases, (jTaf 
and (|7bp. in the f rest frame as 

p f = (m?, 0,0,0), 
p 7 = (£ 7 ,0,0,£ 7 ), 
p T = (E T ,p T sm.9, 0,p T cos 9), 
Plsp = (E,p x ,0,p z ), (8) 

where the z-axis is taken along the photon momentum 
direction, and the j/-axis is along ~p^ 7 x ~p? T , the normal of 
the decay plane. 

Using the generated helicity amplitudes and the above 
kinematical variables, we investigate photon polarizations 
by means of Stokes parameters, P\, P2, and P3, which are 
related with the photon density matrix as 

dp\\> _ 1 / >p \ dr sum 

dE y dcos9 2 I ^ l<7i ' dE y d cos {) 

with the Pauli sigma matrices 0^. dr sum = dp ++ + dp is 

the usual spin-summed differential decay rate. The density 
matrix is calculated as 

d P xx< =^—Y j M x M* xl d^ (10) 

where M.\ is the helicity amplitude with the photon helic- 
ity A, and is the three-body phase space factor. The 
summation symbol implies the summation over the tau 
and gravitino/neutralino helicities. By definition, Stokes 



Diagrams by MadGraph ta1 - -> ta- gro a 




graph 7 



Fig. 1. Feynman diagrams for the radiative stau decay process 
for the G LSP case, f — > tC?7, generated by MadGraph. tal, 
ta, G~, A, and Ni denote a stau, a tau-lepton, a gravitino, a 
photon, and neutralinos, respectively. 

Diagrams by MadGraph ta1 - -> ta- n1 a 
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graph 1 graph 2 

Fig. 2. The same as Fig.[T] but for the Xi LSP case, f — > rxi7- 
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(a) gravitino LSP 

(b) neutralino LSP 
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Fig. 3. Angular dependence of the Stokes parameters of the 
radiated photon for the f decay process, tr — > tGj (a) and 
tr — > Txlj (b), where 9 is the decay angle between the photon 
and the tau-lepton. We set m,f = 150 GeV, mLSP = 75 GeV 
and E-, = 40 GeV. 



parameters take real values from —1 to 1, and P3 shows 
the right-left asymmetry of circular polarizations, while 
Pi and P2 present linear polarizations, which reflect the 
interference between the amplitudes for the right- and left- 
handed photons. 

In Fig. [3l we show the cos 9 dependence of the Stokes 
parameters of the radiated photon for tr — !• rGj (a) and 
fji — > TXi7 (b), where we use 

m? = 150 GeV and m LSP = 75 GeV, (11) 

and fix the photon energy at 



E 1 = 40 GeV. 



(12) 



For the G LSP scenario (a), we take four neutralino masses 
(200, 250, 300, 350) GeV as an example. The 



as m,po 



Pi 



P3 is also shown 



degree of polarization P = \fPy 
with a thick line. Radiated photons are almost fully po- 
larized (P ~ 1) for the both LSP scenarios, except around 
cos6> = —0.95 for the G LSP scenario, where photons are 
close to being unpolarized (P ~ 0) . 

In the cos 9 > region, the photon bremsstrahlung 
amplitude (graph 2 in Figs. [1] and [2]) is dominant and the 
G-LSP and x?-LSP cases are very similar since only ±1/2- 
hclicity states of the gravitino are allowed. In the cos 9 < 
region, on the other hand, the neutralino propagating am- 
plitudes and the four-point interaction amplitude, graph 
3 to 7 in Fig. [TJ become important, which allow the grav- 
itino to take ±3/2 helicities as well. Note that the ampli- 
tude corresponding to the graph 1 in Figs. 1 and 2 always 
vanishes. Since the gravitino has the large mass in this ex- 
ample, spin-3/2 components dominate spin-1/2 ones, and 
P3 for the G LSP shows distinct behavior from those for 
the xi LSP. Especially, for cos 9 ~ —1, the difference is 
significant; P3 = —0.8 (almost left-handed photon) for 
the G LSP, while P3 = +1 (right-handed photon) for the 
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Fig. 4. Missing invariant mass distributions for e + e~ — > 
X1X1 ~~ > 77GG at yfs — 190 GeV. The cases for the neutralino 
mass m x — 75 and 90 GeV are shown as a solid and dashed 
line, respectively, with the normalized cross section after kine- 
matical cuts of (1141) . 



Xi LSP. Those behavior holds for heavier neutralinos and 
agrees with the results of Ref. [TU], where the neutralino 
intermediate diagrams are neglected. 

Since the photon helicity measurements require a po- 
larized detector, we also examine linear polarizations Pi 
and P2. In both scenarios, the linear polarization perpen- 
dicular to the decay plane vanishes (P2 = 0), and Pi tends 
to behave similarly, but slightly larger |Pi| is expected in 
the backward direction (cos 9 < 0) for the gravitino LSP 
case (a). 



2.2 Neutralino NLSP 

As a sample result for the neutralino NLSP scenario, we 
consider the process 



XIX? 



(jG)( 7 G) -> 77#- 



(13) 



Figure [4] shows the distributions of the missing invariant 
mass at y/s = 190 GeV for the neutralino mass m x = 75 
and 90 GeV with the normalized cross section after kine- 
matical cuts. The gravitino mass is fixed at an eV order so 
that xi decays instantly without leaving the production 
point. Here we use the same cuts as in Ref. [TTj : 

I cos# 7 | < 0.95, pt 7 > 0.065 P boam , (14a) 
0.2 < P 7 /P b cam < 0.8, (14b) 



with Pbo 
inHU. 



= s/s/2, and our results agree well with Fig. 16 
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3 Summary 

In this paper, we have added new HELAS subroutines to 
calculate helicity amplitudes with massive spin-3/2 par- 
ticles (massive gravitinos) to the HELAS library. They are 
coded in such a way that arbitrary amplitudes with exter- 
nal gravitinos can be generated automatically by MG, after 
slight modifications. All the codes have been tested care- 
fully by making use of the gauge invariance of the helicity 
amplitudes. 
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Table 1. List of the new vertex subroutines in HELAS system. 



subroutine computes the wavefunction with the w-spinor(3 
The output RI(18) is a complex 18-dimensional array, 
among which the first 16 components contain the wave- 
function as 

RI(4/x+i) = R(// + l,i), (15) 

namely 



A HELAS subroutines for spin-3/2 particles 

In this appendix, we list the contents of all the new HELAS 
subroutines that are needed to evaluate processes based 
on the effective Lagrangian of ([2]) with external spin-3/2 
gravitinos. 

To begin with, in Add. I A. li the subroutines to compute 
external lines for a massive spin-3/2 particle are presented. 
Next, in Apps. IA.2l to [A~5l we explain vertex subroutines 
listed in Table [TJ which compute interactions of a gravitino 
with SM and SUSY particles. Finally, we briefly mention 
how we test our new subroutines in App. IA.6I 



A.l Spin-3/2 wavefunction 

A. 1.1 IRXXXX 

This subroutine computes the flowing-In Rarita-Schwinger 
(R-S) spin-3/2 wavefunction; namely 4>£{p, A) and ^{p 1 A), 
in terms of its four-momentum p and helicity A, and should 
be called as 

CALL IRXXXX(P,RMASS,NHEL,NSR , Rl). 

The input P(0:3) is a real four-dimensional array which 
contains the four-momentum p^ of the spin-3/2 particle, 
RMASS is its mass, NHEL (= ±3,±1) specifies its helicity A 
in unit of 1/2, and NSR specifies whether the fermion is 
particle or anti-particle. If NSR = 1 the fermion is particle 
and the subroutine computes the wavefunction with the 
u-spinor. If NSR = -1 the fermion is anti-particle and the 



RI( 1) = 


R(l,l), 


RI( 2) = 


R(l,2), 


RI( 3) = 


R(l,3), 


RI( 4) = 


R(l,4), 


RI(16) = 


R(4,4), 



where 



R( M +l,i) 



<(P.A) 



for NSR = 1, 
for NSR = -1. 



(16) 



Here, i = 1, 2, 3, 4 denotes each u- or w-spinor component. 
The last two of RI (18) contain the four- momentum along 
the fermion number flow, 



(RI(17), RI(18)) = NSR(P(0) + iP(3), P(l) + iV{2)). 



by 



(17) 



When the four-momentum of the R-S fermion is given 



p^ = (E, \p\ sin6>cos</), |p| sin#sin(/), \p\cos9), (18) 



Although the gravitino is a Majorana particle, the HELAS 
convention requires both of the u- and u-spinors for the calcu- 
lations of amplitudes and their proper interference; see App. A 
in [12]. 
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its helicity states can be expressed as 
$f(p,+3/2) = e"(p,+)u(p,+), 



where 



^(p,+l/2) = A /-^(p,0)«(p 1 +) 



f \l-e> i (p,+)u(p,-)e* 



+ y- e /*(p,Q)u(p,-)e** 
^(p,-3/2) = e»( J >,-)u(p,-)e i *, 



(19) 



by using the vector boson wavefunctions e^(p, A) and the 
spinor wavefunctions u(p, A) that obey the relations 



J^u{p,+) = e i *u{p,-), 



(20) 
(21) 



where J_ = J x — iJ y is the J z lowering operator. The vec- 
tor and spinor wavefunctions in the HELAS convention [3] 
satisfy above relations. Similarly, ip%(p, A) is given by the 
v-spinors and the conjugated vector wavefunctions. The 
above helicity states satisfy the irreducibility conditions 
and the Dirac equation, 

7 ^(P,A)=0, p„C(M)=0, (22) 
(j*-m 3 /2X(p,A)=0, (23) 

and the completeness relation is 

+3/2 

E ^feA)^(p,A) = P^(p), (24) 

A=-3/2 



where 
with 



(p) = {i>+m m )[n^{p)+ l -n^{p)n^(p) lalp 



II»»(p) = -gT + 

TO 3/2 



(25) 
(26) 



A. 1.2 DRXXXX 



This subroutine computes the flowing-Out R-S wavefunc- 
tion; namely, il>u(Pi A) and ip v (p,X), and should be called 
as 

CALL ORXXXX(P,RMASS,NHEL,NSR , RO). 

As in the subroutine IRXXXX, the output R0(18) is a com- 
plex 18-dimensional array, among which the first 16 com- 
ponents contain the wavefunction as 



^ I %.{p,X) for NSR = -1, V ' 



and the last two arc the four-momentum 



(R0(17), R0(18)) = NSR (P(0) + iP(3), P(l) + iP(2)). 



(29) 



A. 2 FRS vertex 

The FRS vertices are obtained from the interaction La- 
grangian among a fermion, a R-S fermion and a scalar 
boson: 



£frs = -ii^ 7 V[GR(l)^i + GR(2)P«]/a,S* + h.c. 



(30) 



with the notation i? p = and the chiral-projection 

operator P R/L = |(1 ±75). GR(1) and GR(2) are the 
relevant left and right coupling constants. For instance, in 
the case of the quark-gravitino-squark interaction, q-G-q a , 
those couplings are 



GR(1) = GFRSL(l) = GFRS 
GR(2) = GFRSL(2) = 

'GR(l) = GFRSR(l) = 
GR(2) = GFRSR(2) = -GFRS 



where 



A. 2.1 I0RSXX 



GFRS = 1/V2M 



PI- 



for a = L, (31) 
for a = R, (32) 

(33) 



This subroutine computes an amplitude of the FRS vertex 
from wavefunctions of a flowing-In fermion, a flowing-Out 
R-S fermion and a Scalar boson, and should be called as 

CALL I0RSXX(FI,R0,SC,GR , VERTEX). 

The input FI (6) is a complex six-dimensional array which 
contains the wavefunction of the flowing-In Fermion and 
its four-momentum as 



P' 



(5ReFl(5), 5ReFl(6), 3mFl(6), 3mFl(5)). 



The input RO ( 18) is a complex 18-dimensional array which 
consists of the wavefunction and the four-momentum of 
the flowing-Out R-S fermion; see the ORXXXX subroutine 
in App. IA.1.21 while the input SC(3) is a complex three- 
dimensional array which contains the wavefunction of the 
Scalar boson, SC(1), and its four-momentum as 



R0(4^i+ i) =R(/x + l,i), 



(27) 



(5fteSC(2), 5fteSC(3), 3mSC(3), 3mSC(2)). 
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The input GR(2) is the complex coupling constant, such 
as in and ([32]) in units of GeV^ 1 . The output VERTEX 
is a complex number in units of GeV: 

VERTEX = (R0) M SC(1)#7 M [GR(1)P L + GR(2)Pr](FI), (34) 

where we use the notations 

fFl(l)\ 

(35) 

VFI(4), 

(R0) M = (RQ(4/x + 1), RD(4/i + 2), R0(4^i + 3), R0(4^, + 4)). 

(36) 

A. 2. 2 IRDSXX 

This subroutine computes an amplitude of the FRS vertex 
from wavefunctions of a flowing-In R-S fcrmion, a flowing- 
Out fermion and a Scalar boson, and should be called as 

CALL IRDSXX (RI,F0,SC,GR , VERTEX). 

The input RI (18) is a complex 18-dimensional array which 
consists of the wavefunction and the four-momentum of 
the flowing-In R-S fermion; see the IRXXXX subroutine 
in App. IA.1.11 while the input F0(6) is a complex six- 
dimensional array which contains the wavefunction of the 
flowing-Out Fermion and its four-momentum as 

p^ = (3?eF0(5),3?eF0(6),3mF0(6),3mF0(5)). 

The output VERTEX is a complex number: 

VERTEX = -(F0)SC(1)[GR(1)*Pr + GR(2)*P L ]7^(Rl) M ,, 

(37) 

where is the momentum of the scalar boson and we use 
the notations 



(38) 



Rl(4/x + l) 
Rl(4/i + 2) 
Rl(4^i + 3) 
Rl(4 M + 4) 

(FD) = (F0(1),F0(2),F0(3),F0(4)). (39) 



and 

FS0R(5) = R0(17) + SC(2), 
FSDR(6) = R0(18) + SC(3). 

Here we use the notation 



(41) 
(42) 



(FSOR) = (FSDR(l), FS0R(2), FSDR(3), FS0R(4)), (43) 

and p is the momentum of the off-shell fermion given in 
and (@5J as 

p M = (3tfeFS0R(5),9tfeFS0R(6),3TOFS0R(6),3TOFS0R(5)). 



A. 2. 4 FSIRXX 

The subroutine computes an off-shell Fermion wavefunc- 
tion made from the interaction of a Scalar boson and a 
flowing-In R-S fermion by the FRS vertex, and should be 
called as 

CALL FSIRXX (RI , SC , GR , FMASS , FWIDTH , FSIR). 
The output FSIR(6) is a complex six-dimensional array: 



(FSIR) = -- 



i{i> + rap) 



m. 



impTp 



SC(1) 



x [iGR(l)*P fi + zGR(2)*P L ] 7 ^(Rl) M , (44) 



and 



FSIR(5) = RI(17) - SC(2), 
FSIR(6) = RI(18) - SC(3). 

Here we use the notation 



(FSIR) = 



and the momentum p is 



/FSIR(l) 
FSIR(2) 
FSIR(3) 

\FSIR(4) 



(45) 
(46) 



(47) 



P 



(3?eFSIR(5), 5ReFSIR(6), 3toFSIR(6), 3toFSIR(5)). 



A.2.3 FSDRXX 



A. 2. 5 HIORXX 



This subroutine computes an off-shell Fermion wavefunc- 
tion made from the interaction of a Scalar boson and a 
flowing-Out R-S fermion by the FRS vertex, and should be 
called as 

CALL FSORXX(RD,SC,GR, FMASS, FWIDTH , FSOR), 

where FMASS and FWIDTH are the mass and the width of 
the fermion, mp and Pp. The output FSOR (6) gives the 
off-shell fermion wavefunction multiplied by the fermion 
propagator and its four-momentum, which is expressed as 
a complex six-dimensional array: 



(FSOR) = (RO)„SC(l)#y"[zGR(l)Pi 



zGR(2)P fl ] 
i + mp) 



impTp ' 



(40) 



This subroutine computes an off-shell scalar current H 
made from the interaction of a flowing-In fermion and 
a flowing-Out R-S fermion by the FRS vertex, and should 
be called as 

CALL HIORXX (FI , R0 , GR , SMASS , SWIDTH , HIOR), 

where SMASS and SWIDTH are the mass and the width of 
the scalar boson, mj and Ps. The output HIOR (3) gives 
the off-shell scalar current multiplied by the scalar boson 
propagator and its four-momentum, which is expressed as 
a complex three-dimensional array: 

HI0R(1) = — 



q 2 - nig + imsTs 
x (R0) M ^y[iGR(l)Pi + iGR(2)P fl ](Fl), (48) 
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and 



A. 3.1 IORVXX 



HI0R(2) = -FI(5) + R0(17), 
HI0R(3) = -FI(6) + R0(18). 



(49) 
(50) 



The momentum q is 

= (3?eHI0R(2),KeHI0R(3),SmHI0R(3),$mHI0R(2)). 

A. 2. 6 HIROXX 

This subroutine computes an off-shell scalar current H 
made from the interaction of a flowing-In R-S fermion and 
a flowing-Out fermion by the FRS vertex, and should be 
called as 

CALL HIROXX (RI , F0 , GR , SMASS , SWIDTH , HIRO). 
The output HIRO (3) is a complex three-dimensional array: 



This subroutine computes an amplitude of the FRV vertex 
from wavefunctions of a flowing-In fermion, a flowing-Out 
R-S fermion and a Vector boson, and should be called as 

CALL IORVXX (FI,R0,VC,GR , VERTEX). 

The input VC(6) is a complex six-dimensional array which 
contains the Vector boson wavefunction and its momen- 
tum as 

<f = (5ReVC(5),3?eVC(6),3mVC(6),3mVC(5)). 

The input GR is the coupling constant in (|55[) . The output 
VERTEX is a complex number: 

VERTEX = (ROy^, f]y[GR(l)P z + GR(2)P S ](FI), (56) 

where we use the notation 



= vc(/i + 1) 



(57) 



HIR0(1) 



q 2 — urig + imsPs 
x (F0)[tGR{i)*P B + iGR(2)*P L ] 7 ^(Rl) M , (51) 



and 



HIR0(2) = -RI(17) + F0(5), 
HIR0(3) = -RI(18) + F0(6). 



(52) 
(53) 



The momentum q is 

= (3f?eHIR0(2),5ReHIR0(3),Q?mHIR0(3),3mHIR0(2)). 

Before turning to the FRV vertex, it should be noticed 
here that the conventional factors of i in the vertices and 
those in the propagators are both included in the off-shell 
wavefunctions, such as (14T))) above, according to the HELAS 
convention. The HELAS amplitude, obtained by the ver- 
tices, such as ([Ml) , gives the contribution to the T matrix 
element without the factor of i. See more details in the 
HELAS manual [3]. 

A.3 FRV vertex 

The FRV vertices are obtained from the interaction La- 
grangian among a fermion, a R-S fermion and a vector 
boson: 

£ F rv = -m^^Yh^i^PL + gr(2)Pr]/ a v v; + h.c. 

(54) 

We note that, although both a gravitino and a gaugino are 
Majorana in most cases, the Hermitian conjugate term is 
necessary for MG; practically, either the first or second term 
is used in calculations of amplitudes. The corresponding 
coupling constant to the effective Lagrangian of ([2]) is 



A. 3. 2 IROVXX 

This subroutine computes an amplitude of the FRV vertex 
from wavefunctions of a flowing-In R-S fermion, a flowing- 
Out fermion and a Vector boson, and should be called as 

CALL IROVXX (RI,F0,VC,GR , VERTEX). 

The output VERTEX is 

VERTEX = -(F0)[GR(1)*P R + GR(2)*P L ] 7 ^[f, rf(RI) M . 

(58) 

A. 3. 3 FVORXX 

This subroutine computes an off-shell Fermion wavefunc- 
tion made from the interaction of a Vector boson and a 
flowing-Out R-S fermion by the FRV vertex, and should be 
called as 

CALL FVORXX (R0 , VC , GR , FMASS , FWIDTH , FVOR) . 

What we compute here is 

(FVOR) = (ROy^, f]^[iGR(l)P L + iGR(2)P R ] 

x „ ^ + WF) _ , (59) 



p 2 — m F + imF^F 



and 



FV0R(5) = R0(17) + VC(5), 
FV0R(6) = R0(18) + VC(6), 



(60) 
(61) 



GR(1) = GR(2) = GFRV = 1/4M P1 . 



(55) 



where we use the notation 

(FVOR) =(FV0R(1),FV0R(2),FV0R(3),FV0R(4)), (62) 
and the momentum p is 
p M = (3f?eFV0R(5),5ReFV0R(6),3mFV0R(6),3TOFV0R(5)). 
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A. 3. 4 FVIRXX 



A. 3. 6 JIROXX 



This subroutine computes an off-shell Fermion wavefunc- 
tion made from the interaction of a Vector boson and a 
flowing-In R-S fermion by the FRV vertex, and should be 
called as 

CALL FVIRXX (RI , VC , GR , FMASS , FWIDTH , FVIR). 

What we compute here is 



(FVIR) 



p 2 — m? F + imFTp 
x {iGR(i)*P R + iGR(2)*P L ]^[f, (63) 



and 



FVIR(5) = RI(17) - VC(5), 
FVIR(6) = RI(18) - VC(6), 

where we use the notation 

(FVIR(l)N 
FVIR(4), 



(64) 
(65) 



(66) 



and the momentum p is 
p 11 = (KeFVIR(5),5ReFVIR(6),3mFVIR(6),3mFVIR(5)). 



This subroutine computes an off-shell vector current J 
made from the interaction of a flowing-In R-S fermion and 
a flowing-Out fermion by the FRV vertex, and should be 
called as 

CALL JIROXX (RI , FO , GR , VMASS , VWIDTH , JIRO). 
The output JIRO (6) is 



JIR0(i/+ 1) 



q 2 — my + imyTy 



-fr 



p v . 



qPq v 



x (F0)[zGR(l)*P fl + zGR(2)*P L ] 7 ^[7 p , ^(Rl) M (71) 
for the massive vector boson, or 

JIR0(^+ 1) = ^ 

q 2 

x (F0)[iGR(l)*P B + zGR(2)*P L ] 7 ^, ^(Rl) M (72) 

for the massless vector boson, and 

JIR0(5) = -RI(17) + F0(5), (73) 
JIR0(6) = -RI(18) + F0(6). (74) 

Here the momentum q is 

g' 1 = (3fJeJIR0(5),3?eJIR0(6),3mJIR0(6),3mJIR0(5)). 



A. 3. 5 JIORXX 

This subroutine computes an off-shell vector current J 
made from the interaction of a flowing-In fermion and 
a flowing-Out R-S fermion by the FRV vertex, and should 
be called as 

CALL JIORXX (FI,R0,GR, VMASS, VWIDTH , JIOR). 

The input VMASS and VWIDTH are the mass and the width 
of the vector boson, my and IV- The output JI0R(6) 
gives the off-shell vector current multiplied by the vec- 
tor boson propagator and its four-momentum, which is 
expressed as a complex six-dimensional array: 



JI0R(z/+ 1) 



q p q u 



q 2 - my + imyTy \ m v 
x (R0y^ 7p ] 7 ^GR(l)P L +2GR(2)P fi ](Fl) (67) 
for the massive vector boson, or 

JI0R(^+ 1) = 

T 

x (R0) Al [^, 7 ,y ]7^[ l GR(l)P L +iGR(2)P B ](FI) (68) 

for the massless vector boson, and 

JI0R(5) = -FI(5) + R0(17), (69) 

JI0R(6) = -FI(6) + R0(18). (70) 

Here, q is the momentum of the off-shell vector boson, 

q^ =(5ReJIOR(5), 5ReJIOR(6), 3mJI0R(6), 3mJI0R(5)). 

Note that we use the unitary gauge for the massive vector 
boson propagator and the Feynman gauge for the massless 
one, according to the HELAS convention [3]. 



A.4 FRVS vertex 

The FRVS vertices are obtained from the interaction La- 
grangian among a fermion, a R-S fermion, a vector boson 
and a scalar boson: 

£frvs = P ai7 1/ 7 ai [GR(1)Pl + GR{2)P R ]fV;S* + h.c. (75) 

The coupling constant GR is the product of the FRS cou- 
pling constant and the gauge coupling constant of the 
involving gauge boson. For instance, in the case of the 
quark-gravitino-gluon-squark interaction, q-G-g-qL-, those 
couplings are 



GR(1) = GFRGSL(l) = GFRSL(l) * GG(l), 



(76) 



1 GR(2) = GFRGSL(2) = GFRSL(2) * GG(2), 

where GFRSL is defined in (f3"Tj) and GG is the strong cou- 
pling constant 

GG(1) = GG(2) = -g a = -G. (77) 

The sign of the coupling constant is fixed by the HELAS 
convention [3]. 

A.4.1 IORVSX 

This subroutine computes an amplitude of the FRVS vertex 
from a flowing-In fermion, a flowing-Out R-S fermion, a 
Vector boson and a Scalar boson, and should be called as 

CALL IORVSX (FI,RO,VC,SC,GR , VERTEX). 
The output VERTEX gives a complex number: 

VERTEX = (R0) M SC(1) )T 7 A1 [GR(l)P L + GR(2)Pr](FI). 

(78) 
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A.4.2 IRDVSX 



A.4.5 JSIORX 



This subroutine computes an amplitude of the FRVS vertex 
from a flowing-In R-S fcrmion, a flowing-Out fermion, a 
Vector boson and a Scalar boson, and should be called as 

CALL IROVSX(RI,FO,VC,SC,GR , VERTEX). 

The output VERTEX gives a complex number: 

VERTEX = (FO) SC(1)[GR(1)*Pr + GR(2)*P L ]7^ f (Rl) M . 

(79) 

A.4.3 FVSORX 

This subroutine computes an off-shell Fermion wavefunc- 
tion made from the interaction of a Vector boson, a Scalar 
boson and a flowing-Out R-S fermion by the FRVS vertex, 
and should be called as 

CALL FVS0RX(R0, VC, SC, GR, FMASS, FWIDTH , FVSOR). 

The output FVSOR is a complex six-dimensional array: 



(FVSOR) = (RO)^SC(l) Y 7 ^GR(1)P L + iGR(2)P R ] 

+ mp) 
p 2 — m 2 F + impFp 

for the first four components of FVS0R(6) , and 



(80) 



FVS0R(5) = R0(17) + VC(5) + SC(2), (81) 
FVS0R(6) = R0(18) + VC(6) + SC(3), (82) 



for the momentum p. 



AAA FVSIRX 

This subroutine computes an off-shell Fermion wavefunc- 
tion made from the interaction of a Vector boson, a Scalar 
boson and a flowing- In R-S fermion by the FRVS vertex, 
and should be called as 

CALL FVSIRX(RI, VC, SC, GR, FMASS, FWIDTH , FVSIR). 

The output FVSIR is a complex six-dimensional array: 

(FVSIR = — ^5 SC(1 

p 2 - m| + im F r F 

x [*GR(l)*P fl + iGR(2)*P L }^ f (RI) M (83) 

for the first four components of FVSIR (6), and 

FVSIR(5) = RI(17) - VC(5) - SC(2), (84) 
FVSIR(6) = RI(18) - VC(6) - SC(3), (85) 



This subroutine computes an off-shell vector current J 
made from the interaction of a Scalar boson, a flowing-In 
fermion and a flowing-Out R-S fermion by the FRVS vertex, 
and should be called as 

CALL JSIORX(FI,RO,SC,GR,VMASS,VWIDTH , JSIOR). 

What we compute here is 



JSI0R(^+ 1) = T —. — 

q z - m v + imyi v 



-9 



v 



x (RO)^SC(l) 7p7 ^GR(l)P L + iGR(2)P fl ](FI) (86) 
for the massive vector boson, or 

JSI0R(^+ 1) = -4 
Q 

x (R0) M SC(1) 7 VNGR(1)Pl + *GR(2)P fl ](FI) (87) 

for the massless vector boson, and 

JSI0R(5) = -FI(5) + R0(17) + SC(2), (88) 
JSI0R(6) = -FI(6) + R0(18) + SC(3), (89) 

for the momentum q. 
A.4.6 JSIROX 

This subroutine computes an off-shell vector current J 
made from the interaction of a Scalar boson, a flowing-In 
R-S fermion and a flowing-Out fermion by the FRVS vertex, 
and should be called as 

CALL JSIROX (RI , FO , SC , GR , VMASS , VWIDTH , JSIRO). 

What we compute here is 



JSIR0(t/+ 1) 



q 2 - my + im v r v 



q p q v 



x (FO) SC(l)[iGR(l)*P fl + zGR(2)*P L ] 7 %(Rl) AJ (90) 
for the massive vector boson, or 

JSIR0(i^+ 1) = ^ 

q 2 

x (FO) SC(1)[»GR(1)*Pr + iGR(2)*P L ]7'y (RI)„ (91) 

for the massless vector boson, and 

JSIR0(5) = -RI(17) + F0(5) + SC(2), (92) 
JSIR0(6) = -RI(18) + F0(6) + SC(3), (93) 



for the momentum p. 



for the momentum q. 
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A. 4. 7 HVIDRX 



A. 5.1 I0RVVX 



This subroutine computes an off-shell scalar current H 
made from the interaction of a Vector boson, a flowing- In 
fermion and a flowing-Out R-S fermion by the FRVS vertex, 
and should be called as 

CALL HVIDRX(FI, R0, VC, GR, SMASS, SWIDTH , HVIOR). 

What we compute here is 

HVIOR(l) = 



and 



q 2 — m s + irasTs 
x (R0) M y 7 ^GR(l)P L + iGR(2)P fl ](FI), (94) 



HVI0R(2) = -FI(5) + R0(17) + VC(5), (95) 
HVI0R(3) = -FI(6) + R0(18) + VC(6), (96) 

for the momentum q. 
A.4.8 HVIROX 

This subroutine computes an off-shell scalar current H 
made from the interaction of a Vector boson, a flowing-In 
R-S fermion and a flowing-Out fermion by the FRVS vertex, 
and should be called as 

CALL HVIR0X(RI, F0, VC, GR, SMASS, SWIDTH , HVIRO). 

What we compute here is 

HVIRO (1) 



and 



q 2 - m 2 s + imsTs 

x (F0)[£GR(1)*P R + iGR(2)*P L ] 7 " J (RI) M , (97) 

HVIR0(2) = -RI(17) + F0(5) + VC(5), (98) 

HVIR0(3) = -RI(18) + F0(6) + VC(6), (99) 



for the momentum q. 



A. 5 FRVV vertex 

The FRVV vertices are obtained from the interaction La- 
grangian among a fermion, a R-S fermion and two vector 
bosons: 

£ffwv = i/ a6 %[ 7 w ,7 p ]7 M [GR(l)^ + GR(2)P*]/ a V$Vf 

+ h.c. (100) 

with the structure constant f abc , which can be handled 
by the MG automatically. The coupling constant GR is the 
product of the FRV coupling constant and the gauge cou- 
pling constant of the involving gauge boson as in the FRVS 
coupling; see ([75]l. 



This subroutine computes an amplitude of the FRVV vertex 
from a flowing-In fermion, a flowing-Out R-S fermion and 
two Vector bosons, and should be called as 

CALL I0RVVX(FI,R0,VA,VB,GR , VERTEX). 

What we compute here is 

VERTEX = (R0) M [p a , f b blGR(l)PL + GR(2)Pr](FI), 

(101) 

(102) 
(103) 

A. 5. 2 IR0VVX 



where we use the notations 



v a ^ = vaO + i), 
v b ^ = vbO + i). 



This subroutine computes an amplitude of the FRVV vertex 
from a flowing-In R-S fermion, a flowing-Out fermion and 
two Vector bosons, and should be called as 

CALL IR0VVX(RI,F0,VA,VB,GR , VERTEX). 

What we compute here is 

VERTEX = (F0)[GR(1)*Pr + QK{2)* P L ]^[f Q , J b ](Rl) M . 

(104) 



A. 5. 3 FVV0RX 

This subroutine computes an off-shell Fermion wavefunc- 
tion made from the interaction of two Vector bosons and 
a flowing-Out R-S fermion by the FRVV vertex, and should 
be called as 

CALL FVV0RX(R0, VA, VB, GR, FMASS, FWIDTH , FVV0R). 

What we compute here is 

(FVV0R) = (ROyr 1 , FV[iGR(l)P L + zGR(2)P fl ] 



p 2 — nip + impFp 



, (105) 



and 



FVV0R(5) = R0(17) + VA(5) + VB(5), (106) 
FVV0R(6) = R0(18) + VA(6) + VB(6). (107) 



A. 5.4 FVVIRX 

This subroutine computes an off-shell Fermion wavefunc- 
tion made from the interaction of two Vector bosons and 
a flowing-In R-S fermion by the FRVV vertex, and should 
be called as 

CALL FVVIRX(RI, VA, VB, GR, FMASS, FWIDTH , FVVIR). 
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What we compute here is 

+ TYl p ) 



(FVVIR) 



p 2 — m F + impTp 



and 



x [zGR(l)*P fl + iGK{2)*P L )^[f a , H(RI) M , (108) 



FVVIR(5) = RI(17) - VA(5) - VB(5), (109) 
FVVIR(6) = RI(18) - VA(6) - VB(6). (110) 



A. 5. 5 JVIORX 

This subroutine computes an off-shell vector current J 
made from the interaction of a Vector boson, a flowing-In 
fermion and a flowing-Out R-S fermion by the FRVV vertex, 
and should be called as 

CALL JVIORX(FI,RO,VC,GR,VMASS,VWIDTH , JVIOR). 

What we compute here is 



JVI0R(//+ 1) 



q 2 — rriy + irnvTv \ ' m 2 



x (R0y 7p , f ]7"[tGR(l)P L + *GR(2)P fl ](FI) (111) 
for the massive vector boson, or 

JVI0R(i/+ 1) = -4 

T 

x (RO^fr", y]7"[»GR(l)P L + *GR(2)P Jl ](FI) (112) 
for the massless vector boson, and 



JVI0R(5) = -FI(5) + R0(17) + VC(5), 
JVI0R(6) = -FI(6) + R0(18) + VC(6). 



(113) 
(114) 



A.5.6 JVIROX 



This subroutine computes an off-shell vector current J 
made from the interaction of a Vector boson, a flowing-In 
R-S fermion and a flowing-Out fermion by the FRVV vertex, 
and should be called as 

CALL JVIROX (RI ,F0 , VC,GR, VMASS , VWIDTH , JVIRO). 

What we compute here is 



JVIR0O+ 1) = 



q 2 — rriy + imyrv 



-fJ 



pu 



qPq u 



ni 



V 



x (F0)[zGR(l)*P fl + iGR(2)*P L ] 7 ^[ 7p , f] (Rl) M (115) 
for the massive vector boson, or 

JVIR0(z/+ 1) = 

x (F0)[»GR(l)*P fl + iGK(2)*P L ]^[Y, f](Rl) M (116) 

for the massless vector boson, and 

JVIR0(5) = -RI(17) + F0(5) + VC(5), (117) 
JVIR0(6) = -RI(18) + F0(6) + VC(6). (118) 



A. 6 Checking for the new HELAS subroutines 

The new HELAS subroutines are tested by using the gauge 
invariance of the helicity amplitudes. In particular, we use 
the following processes; 

qg -> qG for IORSXX, IROSXX, FSORXX, FSIRXX, 
HIORXX, HIROXX, IORVSX, IROVSX, 

gg -> gG for IORVXX, IROVXX, FVORXX, FVIRXX, 
JIORXX, JIROXX, IORVVX, IROVVX, 

qg -> qGg for FVSORX, FVSIRX, JSIORX, JSIROX, 
HVI0RX,HVIR0X, 

gg -> gGg for FVVORX, FVVIRX, JVIORX, JVIROX. 

More explicitly, we express the helicity amplitudes of the 
above processes as 

Mx^ = %{p & , Xg) T^e v {p g , X g ) (119) 

or 

M XsXg = T*» MPg> Afi) e v {p g ,X g ) (120) 

with an external spin-3/2 and a gluon wavefunction. The 
identity for the ST/ (3) gauge invariance 



p gv TK" = 



(121) 



tests all the above subroutines thoroughly. We also test 
the agreement of the helicity-summcd squared amplitudes 
at arbitrary Lorentz frames. 



B Implementation of spin-3/2 gravitinos into 
MadGraph 

In this appendix, we describe how we implement spin-3/2 
gravitinos and their interactions into MG. 

First, using the default mssm model in MG/MEv4 [5], we 
make our new model directory, mssm_gravitino, includ- 
ing a massive gravitino (particles.dat) and its interac- 
tions with SM and SUSY particles (interactions.dat 



3-point 


coupling 


58 






GR 




FRS 


q 


gro 


qi 




GFRSL 






q 


gro 


qr 




GFRSR 




FRV 


g° 


gro 


g 




GFRV 




4-point 


coupling 


5S 






GR 




FRVS 


q 


gro 


g 


qi 


GFRGSL 


= GFRSL*GG 




q 


gro 


g 


qr 


GFRGSR 


= GFRSR*GG 


FRVV 


g° 


gro 


g 


g 


GGORGG 


= GFRV*G 



Table 2. List of the coupling constants for each gravitino ver- 
tex involving SUSY QCD particles. All the particles and the 
coupling constants are written in the MG notation, gro stands 
for a massive gravitino, q represents a light quark, and ql/qr is 
a left/right-handed squark. g and go are a gluon and a gluino, 
respectively. GR is a non-renormalizable coupling constant de- 
fined in each subroutine in App. [X] 
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and couplings . f); we show the coupling constants for 
each gravitino vertex involving SUSY QCD particles in 
Table [2] as examples. Then we add all the new HELAS sub- 
routines for spin-3/2 gravitinos to the HELAS library in MG. 
Since the present MG does not handle spin-3/2 particles, we 
further modify the codes in MG to tell it how to generate 
the FRS, FRV, FRVS and FRVV type of vertices and helicity 
amplitudes, and how to deal with the helicity of external 
spin-3/2 particles. 
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